home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / source / src / arith / irandom.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-11-16  |  1.2 KB  |  65 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  irandom.c
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. /*     irandom.c    RD, 26.03.90 */
  17.  
  18. #include <LEDA/impl/iint.h>
  19. #include <LEDA/impl/iloc.h>
  20.  
  21. #include <time.h>
  22. #include <sys/types.h>
  23.  
  24. static PLACE Prandom()
  25.     /* return  random 16/31 bits; */
  26. {    unsigned int x;
  27.         time_t seed;
  28.     static BOOLEAN init=FALSE;
  29.  
  30.     if (!init) {
  31.         init=TRUE;
  32.         time(&seed);
  33.         srandom((int)seed);
  34.     }
  35.     return random();
  36. }        /* Prandom */
  37.  
  38. void IasrandomI(a, b)
  39.     Integer *a;
  40.     const Integer *b;
  41. /* waehle a zufaellig mit 0<=|a|<|b|, a->sign=b->sign.
  42.    a und b muessen verschieden sein. */
  43. {    register int i;
  44.     register pPLACE pa;
  45.     register int nl;
  46.     nl=b->length;
  47.     if (nl>a->maxlength) {
  48.         delvec(a->vec, a->maxlength);
  49.         a->maxlength=nl;
  50.         a->vec=newvec(&a->maxlength);
  51.     }
  52.     pa=a->vec;
  53.     pa[nl-1]=Prandom()%(b->vec[nl-1]);
  54.     for (i=nl-2; i>=0; i--)
  55.         pa[i]=Prandom();
  56.     i=nl;
  57.     while (!pa[i-1] && i>0)
  58.         i--;
  59.     a->length=i;
  60.     if (i)
  61.         a->sign=b->sign;
  62.     else
  63.         a->sign=PLUS;
  64. }            /* IasrandomI */
  65.